TABLE 2 

64 Psc Dataset, Best-Fit Model Predictions, and Residuals 



MJD 


Cal Date/UT 


Meas V"^ 


Error 


Fred 


Resid 


WL (^ni) 


u (m) 


V (m) 


HA (hrs) 


50690.3719 


8/30/97 8:55:32 


0.816 


0.031 


0.830 


-0.013 


2.18 


-52.59 


-94.30 


-1.11 


50690.4278 


8/30/97 10:15:59 


0.621 


0.034 


0.661 


-0.040 


2.18 


-33.35 


-98.75 


0.24 


50690.4747 


8/30/97 11:23:34 


0.331 


0.026 


0.331 


0.000 


2.17 


-13.90 


-100.80 


1.37 


50691.3888 


8/31/97 9:19:53 


0.613 


0.082 


0.665 


-0.052 


2.18 


-46.45 


-96.08 


-0.63 


50691.4388 


O /01 /AT 1A 01 C C 

8/31/97 10:31:55 


A O 1 O 
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A AO 1 
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A O O/^ 
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O 1 o 
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AA C O 
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A CT T 

0.57 


50692.3162 


9/1/97 7:35:21 


0.239 


0.019 


0.240 


-0.002 


2.19 


-64.43 


-88.86 


-2.32 


50692.3666 


9/1/97 8:47:56 


0.120 


0.008 


0.121 


-0.001 


2.19 


-52.54 


-94.32 


-1.10 


50692.4367 


a/i /at 1A OO /l/^ 

9/1/97 10:28:46 


0.099 


0.011 


0.097 


0.001 


2.19 


-27.64 


-99.55 


0.58 
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O.OOd 
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0.105 
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07 AO 
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0.62 


50693.2995 


9/2/97 7:11:18 


0.421 


0.018 


0.437 


-0.015 


2.19 


-66.63 


-87.18 


-2.65 


50693.3080 


9/2/97 7:23:27 


0.454 


0.024 


0.452 


0.001 


2.19 


-65.36 


-88.20 


-2.45 


C A/? AO O O 1 A 


A /O /AT T.C^J.OC 


r\ con 
0.530 


0.029 


0.526 


A A A C 

0.005 


O 1 A 

2.19 
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-60.99 


AA O O 
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1 A A 

-1.90 


50693.3705 


A /a /at o r o at 

9/2/97 8:53:27 


0.751 


0.038 


0.727 


0.025 


2.19 


-50.57 


-94.94 


-0.94 


50693.4116 


A /a /at a r" a /i o 

9/2/97 9:52:42 


0.870 


0.059 


0.928 


-0.058 


2.19 


-36.40 


-98.24 


0.05 


50720.3756 


9/29/97 9:00:53 


0.447 


0.039 


0.492 


-0.045 


2.19 


-21.25 


-100.24 


0.95 


C ATO A /I 1 A C 


A /OA /AT A. CI .ACT 

y/zy/y y:oi:Uo 


n o on 
0.839 


A An c 

0.095 


A O /I C 

0.845 


A AA^ 
-0.00 


O 1 A 

2.19 


-D.l ( 


1 A1 1 O 

-101.12 


1. ()) 


50720.4355 


9/29/97 10:27:00 


0.991 


0.128 


0.956 


0.035 


2.19 


4.89 


-101.15 


2.40 


50723.3560 


10/2/97 8:32:39 


0.257 


0.040 


0.209 


0.048 


2.20 


-26.00 


-99.75 


0.68 


50725.2513 


10/4/97 6:01:53 


0.678 


0.044 


0.689 


-0.011 


2.20 


-59.23 


-91.73 


-1.71 


50725.3018 


10/4/97 7:14:38 


0.678 


0.029 


0.691 


-0.013 


2.19 


-44.48 


-96.57 


-0.49 


50725.3579 


10/4/97 8:35:21 


0.780 


0.045 


0.809 


-0.029 


2.19 


-22.96 


-100.08 


0.86 


50726.2670 


10/5/97 6:24:31 


0.797 


0.157 


0.872 


-0.075 


2.20 


-54.47 


-93.65 


-1.27 


50726.3002 


10/5/97 7:12:16 


0.830 


0.155 


0.932 


-0.102 


2.20 


-44.11 


-96.66 


-0.47 


50726.3251 


10/5/97 7:48:12 


0.973 


0.145 


0.955 


0.018 


2.20 


-35.02 


-98.48 


0.13 


50726.3555 


10/5/97 8:31:54 


1.038 


0.228 


0.919 


0.119 


2.20 


-22.83 


-100.09 


0.86 


50753.1680 


11/1/97 4:01:58 


0.582 


0.037 


0.584 


-0.002 


2.19 


-60.80 


-90.97 


-1.87 


50753.1769 


11/1/97 4:14:40 


0.566 


0.034 


0.579 


-0.013 


2.19 


-58.75 


-91.94 


-1.66 


50753.2000 


11/1/97 4:47:56 


0.589 


0.029 


0.589 


0.000 


2.19 


-52.57 


-94.31 


-1.11 


50753.2236 
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11/1/97 5:21:59 
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0.037 


0.628 


-0.006 


2.19 
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-96.43 


-0.54 


50753.2490 


11/1/97 5:58:32 


0.672 


0.035 


0.696 


-0.023 


2.19 


-35.95 


-98.32 


0.07 


50755.1742 


11/3/97 4:10:54 


0.190 


0.016 


0.190 


0.000 


2.19 


-58.06 


-92.25 


-1.59 


50755.1877 


11/3/97 4:30:14 


0.233 


0.022 


0.250 


-0.018 


2.19 


-54.52 


-93.63 


-1.27 


50755.1962 


11/3/97 4:42:29 


0.277 


0.026 


0.293 


-0.016 


2.19 


-52.07 


-94.47 


-1.07 


50755.3236 


-i-t i i ^ A r' ^ A 

11/3/97 7:45:54 


0.940 


0.075 


0.965 


-0.025 


2.19 


-2.36 


-101.19 


2.00 


50755.3314 


11/3/97 7:57:11 


0.942 


0.098 


0.972 


-0.030 


2.19 


1.10 


-101.20 


2.19 


50771.2294 


11 /ia/at r" OA oi 

11/19/97 5:30:21 


0.430 


0.202 


0.494 


-0.064 


2.19 


-24.16 


-99.95 


0.79 


50771.3006 


11/19/97 7:12:53 


0.504 


0.796 


0.628 


-0.124 


2.19 


6.80 


-101.11 


2.50 


50771.3036 


11/19/97 7:17:11 


2.152 


2.081 


0.642 


1.510 


2.19 


8.11 


-101.07 


2.57 


OU / / D.104U 
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U.Uo / 


z.zu 


/I7 'i'^ 


-yo.oo 


n 7n 

-U. ( u 


50776.2221 


11/24/97 5:19:50 


1.000 


0.333 


0.956 


0.044 


2.20 


-21.51 


-100.22 


0.94 


51053.3629 


8/28/98 8:42:30 


0.726 


0.041 


0.849 


-0.123 


2.21 


-56.79 


-92.77 


-1.47 


51053.3809 


8/28/98 9:08:29 


0.741 


0.053 


0.773 


-0.032 


2.21 


-51.75 


-94.57 


-1.04 


51053.4042 


8/28/98 9:42:04 


0.641 


0.033 


0.626 


0.015 


2.21 


-44.25 


-96.63 


-0.48 


51053.4357 


8/28/98 10:27:26 


0.365 


0.023 


0.364 


0.000 


2.21 


-32.64 


-98.86 


0.28 


51053.4546 


8/28/98 10:54:35 


0.202 


0.012 


0.206 


-0.003 


2.20 


-25.05 


-99.86 


0.73 


51056.3316 


8/31/98 7:57:29 


0.769 


0.061 


0.901 


-0.131 


2.20 


-62.16 


-90.25 


-2.03 


51056.3538 


8/31/98 8:29:30 


0.886 


0.072 


0.928 


-0.042 


2.20 


-57.01 


-92.69 


-1.49 


51056.3767 


8/31/98 9:02:27 


0.853 


0.039 


0.943 


-0.090 


2.20 


-50.54 


-94.95 


-0.94 
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TABLE 2 Confiniicd 



MJD 


Cal Date/UT 


Meas 


Error 


Pred y2 


Resid 


WL (/xm) 


u (m) 


V (m) 


HA (hrs) 


51056.3883 


8/31/98 9:19:06 


0.980 


0.043 


0.948 


0.032 


2.20 


-46.86 


-95.98 


-0.66 


51056.4151 


8/31/98 9:57:47 


0.923 


0.047 


0.953 


-0.030 


2.20 


-37.39 


-98.06 


-0.02 


51056.4475 


8/31/98 10:44:25 


0.915 


0.046 


0.954 


-0.040 


2.20 


-24.58 


-99.91 


0.76 


51056.4713 


8/31/98 11:18:41 


0.963 


0.056 


0.953 


0.010 


2.20 


-14.49 


-100.77 


1.33 


51056.4943 


8/31/98 11:51:48 


1.017 


0.057 


0.947 


0.070 


2.20 


-4.43 


-101.17 


1.89 


51078.3139 


9/22/98 7: 


32:01 


0.672 


0.042 


0.719 


-0.047 


2.21 


-51.37 


-94.69 


-1.01 


51078.3431 


9/22/98 8: 


14:03 


0.555 


0.049 


0.595 


-0.040 


2.21 


-41.76 


-97.19 


-0.31 


51078.3702 


9/22/98 8: 


53:05 


0.383 


0.035 


0.423 


-0.040 


2.20 


-31.57 


-99.02 


0.35 


51098.2705 


10/12/98 6 


:29:31 


0.798 


0.041 


0.839 


-0.041 


2.21 


-47.87 


-95.71 


-0.74 


51098.3024 


10/12/98 7 


:15:27 


0.840 


0.036 


0.849 


-0.009 


2.21 


-36.66 


-98.20 


0.03 


51098.3371 


10/12/98 8 


:05:21 


0.861 


0.034 


0.883 


-0.023 


2.21 


-22.83 


-100.09 


0.86 


51098.3611 


10/12/98 8 


:39:56 


0.885 


0.042 


0.918 


-0.033 


2.21 


-12.57 


-100.88 


1.44 


51098.4005 


10/12/98 9 


:36:41 


0.942 


0.042 


0.956 


-0.015 


2.21 


4.78 


-101.16 


2.39 


51099.2762 


10/13/98 6 


:37:45 


0.761 


0.030 


0.785 


-0.024 


2.20 


-45.07 


-96.43 


-0.54 


51099.3223 


10/13/98 7 


:44:04 


0.907 


0.034 


0.923 


-0.016 


2.21 


-27.80 


-99.53 


0.57 


51121.3010 


11/4/98 7:13:23 


0.813 


0.060 


0.777 


0.036 


2.23 


-11.40 


-100.94 


1.51 


51121.3174 


11/4/98 7:37:02 


0.970 


0.093 


0.899 


0.072 


2.23 


-4.18 


-101.17 


1.90 


51121.3203 


11/4/98 7:41:17 


0.930 


0.115 


0.916 


0.014 


2.23 


-2.88 


-101.19 


1.97 


51122.2230 


11/5/98 5:21:04 


0.859 


0.126 


0.931 


-0.072 


2.21 


-41.76 


-97.19 


-0.31 


51122.2685 


11/5/98 6:26:38 


0.579 


0.059 


0.652 


-0.073 


2.21 


-24.09 


-99.96 


0.79 


51122.2827 


11/5/98 6:47:03 


0.466 


0.058 


0.508 


-0.042 


2.22 


-18.12 


-100.51 


1.13 


51122.2850 


11/5/98 6:50:23 


0.620 


0.125 


0.485 


0.134 


2.21 


-17.13 


-100.59 


1.19 


51122.2954 


11/5/98 7:05:22 


0.547 


0.081 


0.376 


0.171 


2.21 


-12.63 


-100.87 


1.44 


51122.3075 


11/5/98 7:22:46 


0.316 


0.031 


0.252 


0.064 


2.21 


-7.35 


-101.09 


1.73 


51127.1510 


11/10/98 3:37:22 


0.881 


0.032 


0.852 


0.029 


2.22 


-59.24 


-91.72 


-1.71 


51127.1655 


11/10/98 3:58:14 


0.853 


0.087 


0.869 


-0.016 


2.21 


-55.56 


-93.25 


-1.36 


51127.1788 


11/10/98 4:17:26 


0.880 


0.048 


0.895 


-0.016 


2.22 


-51.77 


-94.56 


-1.04 


51127.1983 


11/10/98 4:45:35 


0.956 


0.070 


0.927 


0.029 


2.21 


-45.56 


-96.31 


-0.57 


51127.2479 


11/10/98 5:56:55 


0.889 


0.038 


0.941 


-0.053 


2.21 


-26.97 


-99.64 


0.62 


51127.2673 


11/10/98 6:24:53 


0.835 


0.029 


0.887 


-0.052 


2.21 


-18.86 


-100.46 


1.09 


51127.2880 


11/10/98 6:54:41 


0.752 


0.041 


0.776 


-0.024 


2.22 


-9.92 


-101.00 


1.59 


51127.3079 


11/10/98 7:23:25 


0.552 


0.031 


0.631 


-0.078 


2.22 


-1.12 


-101.20 


2.07 


51127.3228 


11/10/98 7:44:49 


0.536 


0.101 


0.496 


0.040 


2.22 


5.44 


-101.15 


2.43 


51133.1696 


11/16/98 4:04:16 


0.767 


0.046 


0.744 


0.022 


2.23 


-49.56 


-95.24 


-0.87 


51133.1845 


11/16/98 4:25:43 


0.690 


0.067 


0.690 


0.000 


2.23 


-44.68 


-96.53 


-0.51 


51133.1995 


11/16/98 4:47:17 


0.635 


0.033 


0.626 


0.009 


2.23 


-39.39 


-97.68 


-0.15 


51133.2106 


11/16/98 5:03:12 


0.569 


0.029 


0.557 


0.012 


2.23 


-35.25 


-98.44 


0.12 


51133.2222 


11/16/98 5:19:58 


0.516 


0.019 


0.479 


0.038 


2.23 


-30.71 


-99.15 


0.40 


51133.2330 


11/16/98 5:35:34 


0.403 


0.023 


0.401 


0.002 


2.23 


-26.33 


-99.71 


0.66 


51133.2449 


11/16/98 5:52:38 


0.331 


0.036 


0.309 


0.022 


2.23 


-21.41 


-100.23 


0.94 
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ABSTRACT 



We report on the determination of the visual orbit of the double-hned spectroscopic 
binary system 64 Piscum with data obtained by the Palomar Testbed Interferometer 
in 1997 and 1998. 64 Psc is a nearly equal-mass double-lined binary system whose 
spectroscopic orbit is well known. We have estimated the visual orbit of 64 Psc from 
our interferometric visibility data. Our 64 Psc orbit is in good agreement with the 
spectroscopic results, and the physical parameters implied by a combined fit to our 
interferometric visibility data and radial velocity data of Duquennoy and Mayor result 
in precise component masses that agree well with their spectral type identifications. In 
particular, the orbital parallax of the system is determined to be 43.29 ± 0.46 mas, 
and masses of the two components are determined to be 1.223 it 0.021 Mq and 1.170 
± 0.018 M0, respectively. 

Nadal et al. put forward arguments of temporal variability in some of the orbital 
elements of 64 Psc, presumably explained by an undetected component in the 
system. While our visibility data does not favor the Nadal temporal variability 
inference, neither is it definitive in excluding it. Consequently we have performed both 
high dynamic-range near-infrared imaging and spectroscopy of potential additional 
companions to the 64 Psc system. Our imaging and spectroscopic data do not support 
the conjecture of an additional component to 64 Psc, but we did identify a faint object 
with unusual red colors and spectra. 
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1. Introduction 

64 Piscum (HD 4676) is a nearby, short-period (13.8 d) binary system with nearly-equal mass 
{q ~ 0.96) components of spectral type F8 V. Further, the 64 Psc binary is the bright component 



in an apparent visual triple system (Washington Double Star Catalog ~ WDS, Worley & Douglass 



1997 ). The WDS lists two additional visual components B and C with V of 12.6 and 13.0 at 77 



and 71 arcsec separations from 64 Psc, respectively. If these objects are in fact bound to 64 Psc, 



with a system distance of approximately 24 pc as determined by Hipparcos (ESA 1997, Perryman 
|et al. 19"97| ), the minimum physical separations of the dim B and C components to 64 Psc are 
approximately 1800 AU, indicating orbital periods of ~ 5 x 10^ yr. 

Abt and Levy (1976, herein AL76) determined the first double-lined spectroscopic orbit for 
the 64 Psc system. Both Nadal et al (1979, herein N79) and later Duquennoy and Mayor (data 
in 1991a - herein DM91a, and solution in 1991b - herein DM91b) have determined additional 
double-lined orbital solutions for the 64 Psc system, which are all in reasonable agreement with 
each other. Particularly notable is Nadal's argument for the potential presence of an additional, 
undetected object in the 64 Psc system from possible time variation in the system velocity, the 
time of periastron passage, and the argument of periastron. N79 used 1918 data from Lick and 
Mt. Wilson, reanalyzed the AL76 observations, and presented their own data in support of the 
variation hypothesis. Presumably these time-variation effects could not be due to either the B or 
C visual companions listed in the WDS; N79 estimates the period of these effects to be ~ 2700 yr. 

Herein we report the determination of the 64 Psc visual orbit from near-infrared, long-baseline 
interferometric measurements taken with the Palomar Testbed Interferometer. PTI is a 110-m 
K-hand (2.2 ^um) interferometer located at Palomar Observatory, and described in detail elsewhere 
( Colavita et al. 1999a| ). PTI has a minimum fringe spacing of roughly 4 mas at the sky position of 



64 Psc, making this binary system readily resolvable. 
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2. Observations 

The interferometric observable used for these measurements is the fringe contrast or visibility 
(squared) of an observed brightness distribution on the sky. Normahzed in the interval [0:1], a 
single star exhibits visibility modulus in a uniform disk model given by: 

2J,{7tB9/X) 

ttBO/X ^ ' 

where Ji is the first-order Bessel function, B is the projected baseline vector magnitude at the 
star position, is the apparent angular diameter of the star, and A is the center-band wavelength 
of the interferometric observation. The expected squared visibility in a narrow pass-band for a 
binary star such as 64 Psc is given by: 



+ Vi r'^ + 2ViV2r cos(^ B • s) 
(l + r)2 



where Vi and V2 are the visibility moduli for the two stars separately as given by Eq. |T], r is the 
apparent brightness ratio between the primary and companion, B is the projected baseline vector 
at the system sky position, and s is the primary-secondary angular separation vector on the plane 
of the sky ( Pan et al. 199C , Hummel et al. 19951 ). The observables used in our 64 Psc study are 



both narrow-band V'^ from individual spectral channels (Colavita et al. 1999a), and a synthetic 
wide-band V'^, given by an incoherent SNR- weighted average of the narrow-band channels in 
the PTI spectrometer ( polavita 1999b ). In this model the expected wide-band observable is 



approximately given by an average of the narrow-band formula over the finite pass-band of the 
spectrometer: 

1 " 

^J6 = -EKV.(A0 (3) 

where the sum runs over the channels covering the infrared K-band (2 - 2.4 /im) of the PTI 
spectrometer ([Colavita et al. 1999a ). Separate calibrations and hypothesis fits to the narrow-band 



and synthetic wide-band V'^ datasets yield statistically consistent results, with the synthetic 
wide-band data exhibiting superior fit performance. Consequently we will present only the results 
from the synthetic wide-band data. 
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Object 


Spectral 


Star 


64 Psc 


Adopted Model 


Name 


Type 


Magnitude 


Separation 


Diameter (mas) 


HD 166 


KO V 


6.1 V/4.1 K 


16° 


0.66 ± 0.06 


HD 3651 


KO V 


5.9 V/3.9 K 


4.9° 


0.79 ± 0.05 


HD 4628 


K2 V 


5.7 V/3.5 K 


12° 


0.92 ± 0.05 



Table 1: PTI 64 Psc Calibration Objects Considered in our Analysis. The relevant parameters for 
our two calibration objects are summarized. The apparent diameter values are determined from 
effective temperature and bolometric flux estimates based on archival broad-band photometry, and 
visibility measurements with PTI. 



64 Psc was observed in conjunction with our selected calibrator list by PTI on 21 nights 
between 30 Aug 1997 and 16 Nov 1998, covering roughly 32 periods of the system. 64 Psc, 
along with calibration objects, was observed multiple times during each of these nights, and 
each observation, or scan, was approximately 120 sec long. For each scan we computed a mean 
V"^ value from the scan data, and the error in the V'^ estimate from the RMS internal scatter 
( Colavita 1999b| ). 64 Psc was always observed in combination with one or more calibration 
sources within ~ 20 degrees on the sky. For our study we have used three late-type main sequence 
stars as calibration objects: HD 166, HD 3651, and HD 4628. Table || lists the relevant physical 
parameters for the calibration objects. 

The calibration of 64 Psc V'^ data is performed by estimating the interferometer system 
visibility (Vgyg) using calibration sources with model angular diameters, and then normalizing the 
raw 64 Psc visibility by V^y^ to estimate the V'^ measured by an ideal interferometer at that epoch 
(Mozurkewich et al. 1991, Boden et al. 1998| ). Calibrating our 64 Psc dataset with respect to the 
three calibration objects listed in Table || results in a total of 88 calibrated scans on 64 Psc over 
21 nights in 1997 and 1998. 
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3. Orbit Determination 

The estimation of the 64 Psc visual orbit is made by fitting a Keplarian orbit model with 
visibilities predicted by Eqs. || and ^ directly to the calibrated (narrow-band and synthetic wide- 
band) V'^ data on 64 Psc (see Armstrong et al. 1992b, Hummel 1993, Hummel et al. 1995| , Boden 



let al. 1999D . The fit is non-linear in the Keplarian orbital elements, and is therefore performed by 



non- linear least-squares methods (i.e. the Marquardt-Levenberg method. Press et al. 1992 ). As 
such, this fitting procedure takes an initial estimate of the orbital elements and other parameters 
(e.g. component angular diameters, brightness ratio), and evolves that model into a new parameter 
set which best fits the data. However, the chi-squared surface has many local minima in addition 
to the global minimum corresponding to the true orbit. Because the Marquardt-Levenberg method 
strictly follows a downhill path in the manifold, it is necessary to thoroughly survey the space 
of possible binary parameters to distinguish between local minima and the true global minimum. 
In addition, as the V'^ observable for the binary (Eqs. |2| and |^ is invariant under a rotation of 
180°, we cannot differentiate between an apparent primary/secondary relative orientation and 
its mirror image on the sky. Consequently there remains a 180° ambiguity in our determination 
of the longitude of the ascending node, 0, which we quote by convention in the interval [0:180). 
By similar arguments our V'^ observable does not distinguish the longitude of periastron (a;) for 
the primary and secondary component. We have quoted our estimate constrained to be grossly 
(within 180°) consistent with the consensus LUi value of approximately 200° (AL76, N79, DM91b). 

In addition to our PTI visibility data we have used the double-lined radial velocity data from 
DM91a. To incorporate this data we have upgraded our orbit estimation software to utilize both 
interferometric visibility and radial velocity data either separately or simultaneously (see similar 
remarks in Hummel et al. 1998| ). The w-degeneracy discussed above is resolved by the inclusion 



of radial velocity data in our orbital solution (Table |^ , however the determination of Q remains 
ambiguous by 180°. 

In the case of 64 Psc the parameter space is significantly narrowed by the high-quality 
spectroscopic orbit from both N79 and DM91b, and the Hipparcos distance determination sets 



-7- 



the rough scale of the semi-major axis (ESA 1997). Further, because the 64 Psc system is nearly 



equal-mass {q = 0.959 it 0.004 - DM91b), and because interferometric V'^ measurements are 
only weakly sensitive to variations in intensity ratio near equal intensity, we have constrained the 
intensity ratios in our fits to 0.90 [AK = 0.11) as indicated by our component mass solutions (§^, 
and the empirical mass-luminosity relationship of Henry and McCarthy (1992). Finally, at the 
distance of 64 Psc the apparent diameters of the two main-sequence components of the system are 
not strongly resolved by PTI, so we have constrained the estimated diameters of both components, 
to 0.5 mas for the primary component and 0.45 mas for the secondary component, as indicated 
by our system distance and component mass estimates (^Q), and a fit to empirical mass-radius 
data for main sequence eclipsing systems given in Andersen (1991). Given this limited parameter 
space, the correct orbit solution is readily obtained by exhaustive search. 

Figure [l| depicts the relative visual orbit of the 64 Psc system, with the primary component 
rendered at the origin, and the secondary component rendered at periastron. We have indicated 
the phase coverage of our V'^ data on the relative orbit with heavy lines; our data samples most 
phases of the orbit well, leading to a reliable orbit determination. 

Table ^ lists the complete set of measurements in our 64 Psc dataset and the prediction 
based on the best-fit orbit model (our full-fit model. Table ^ for 64 Psc. Table ^ gives V'^ 
measurements and times, measurement errors, model predictions, the photon-weighted average 
wavelength, u — v coordinates, and on-target hour angle for each of our calibrated 64 Psc 
observations. Figure ^a shows four consecutive nights of PTI V"^ data on 64 Psc (30 Aug - 2 Sept 
1997), and V"^ predictions based on the best-fit model for the system (our full- fit model. Table ^). 
The model predictions are in good agreement with the observed data, with an RMS V'^ residual 
of 0.034 (average absolute V'^ deviation of 0.053), and a per Degree of Freedom (DOF) of 0.85. 
Figure |2|b gives a radial velocity phase plot of the DM91a radial velocity data and the predictions 
of our full-fit orbital solution. A radial velocity fit residual histogram is also included. 

Spectroscopic orbit parameters (from N79 and DM91b) and our visual and spectroscopic 
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Fig. 1. — Visual Orbit of 64 Psc. The relative visual orbit of 64 Psc is shown, with the primary and 
secondary objects rendered at Tq (periastron). The heavy lines along the relative orbit indicate 
areas where we have orbital phase coverage in our PTI data (they are not separation vector 
estimates); our data samples most phases of the orbit well, leading to a reliable orbit determination. 
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Fig. 2. — V'^ and RV Fit of 64 Psc. a: Four consecutive nights (30 Aug - 2 Sept 1997) of calibrated 
V'^ data on 64 Psc, and V'^ predictions from the best-fit model for the system. The lower frame 
shows individual V'^ residuals between the calibrated data and best-fit model, b: A phase plot 
of radial velocity data from DM91a and fit predictions from our Full-Fit solution (Table ^). Also 
included is a histogram of RV residuals around the double-lined solution. 



- 10 - 



orbit parameters of the 64 Psc system are summarized in Table |3[ We give the results of separate 
fits to only our V'^ data (our V^-only Fit solution), and a simultaneous fit to our data 
and the double-lined radial velocities from DM91a (our Full-Fit solution) - both with relative 
intensity and component diameters constrained as noted above. For the orbit parameters we 
have estimated from our visibility data we list a total one-sigma error in the parameter estimate, 
and the separate one-sigma errors in the parameter estimates from statistical (measurement 
uncertainty) and systematic error sources. In our analysis the dominant forms of systematic 
error are: (1) uncertainties in the calibrator angular diameters (Table |l|); (2) uncertainty in our 
center-band operating wavelength (Aq ~ 2.2 /im), which we have taken to be 20 nm (~1%); (3) 
the geometrical uncertainty in our interferometric baseline ( < 0.01%); and (4) uncertainties in 
orbital parameters we have constrained in our fitting procedure (i.e. the component intensity ratio 
and angular diameters in both solutions). Different parameters are affected differently by these 
error sources; e.g. our estimated uncertainty in the 64 Psc orbital inclination is dominated by 
measurement uncertainty, while the uncertainty in the angular semi-major axis is dominated by 
uncertainty in the wavelength scale. Conversely, we have assumed that all of the error quoted by 
N79 and DM91b in the 64 Psc spectroscopic parameters is statistical. 

There is a significant (~ 4.2-sigma) difference in the period as estimated by DM91b and 
the period estimates of N79 and PTI. For the full-fit solution our period estimate exploits the 
roughly 7300-day (530 64 Psc orbital cycles) time baseline between the DM91a RV data and our 
own visibility data. Both our period estimates agree well with the N79 determination (the full-fit 
solution agreement with N79 is ^ 0.9-sigma). 

4. Physical Parameters 

Physical parameters derived from our 64 Psc full- fit visual/spectroscopic orbit are summarized 
in Table |^. As in Table ^, for physical parameters we have estimated we quote total sigma 
error, and its statistical and systematic contributions. Notable among these is the high-precision 
determination of the component masses for the system, a virtue of the precision of the DM91a 
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KJl. Ultdi 

Parameter 


N7Q 




PTI 97/98 


V -only r it 


r uli r It 


Feriod (dj 


13.824581 


13.8318 


13.82407 


13.824621 




± 4.0 X 10~^ 


it 1.7 X 10~^ 


zb 4.4 (3.8/2.2) X 10~^ 


zb 1.7 (1.7/0.4) X 10"^ 


To (MJD) 


41933.702 


43468.308 


50905.934 


50905.984 




lb 0.074 


± 0.018 


lb 0.068 (0.048/0.048) 


zb 0.015 (0.015/0.003) 


e 


0.243 ± 0.010 


0.238 ± 0.002 


0.2348 


0.2376 








zb 0.0052 (0.0035/0.0039) 


zb 0.0012 (0.0011/0.0005) 


Ki (km s ^) 


57.53 ± 0.74 


57.31 lb 0.19 




57.35 zb 0.31 (0.31/0.01) 


K2 (km s~^) 


58.77 ± 0.75 


59.77 lb 0.18 




59.95 zb 0.32 (0.31/0.02) 


wi (deg) 


199.6 ± 2.2 


202.3 lb 0.5 


202.2 =b 1.9 (1.3/1.3) 


203.56 zb 0.35 (0.33/0.10) 


Oi (deg) 






63.43 zb 0.68 (0.66/0.16) 


63.60 zb 0.82 (0.81/0.16) 


i (deg) 






74.37 zb 0.96 (0.92/0.27) 


73.80 zb 0.92 (0.91/0.17) 


a (mas) 






6.516 ± 0.061 (0.016/0.059) 


6.527 ± 0.061 (0.015/0.059) 


AiC (mag) 






0.11 


0.11 








0.85 


1.2 (0.95 FV2.8 RV) 








0.053 


0.053 


\RRy \ (km s-i) 




0.65 




0.65 



Table 3: Orbital Parameters for 64 Psc. Summarized here are the apparent orbital parameters 
for the 64 Psc system as determined by N79, DM91b, and PTI. We give two separate fits to our 
data, with and without including the double-lined DM91a radial velocities in the fit (constrained 
parameters are listed in italics). For parameters we have estimated from including our PTI 
observations we separately quote one-sigma errors from both statistical and systematic sources, 
given as {ustat/ <^sys)i and the total error as the sum of the two in quadrature. We have quoted 
the longitude of the ascending node parameter (O) as the angle between local East and the orbital 
line of nodes measured positive in the direction of local North. Due to the degeneracy in our 
observable there is a 180° ambiguity in 17; by convention we quote it in the interval of [0:180). We 
quote mean absolute V'^ and RV residuals in the fits, \RvA ^^'^ respectively. 
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radial velocity measurements on both components and the high inclination of the orbit. We 
estimate the masses of the F8V primary and secondary components as 1.223 ± 0.021 Mq and 
1.170 lb 0.018 Mq, respectively. Remarkably, these are very near the values estimated by N79, 
based on their radial velocities and the relative magnitude in the system. 

The Hipparcos catalog hsts the parallax of 64 Psc as 41.80 ± 0.75 mas ( |ESA 1997] ). The 
distance determination to 64 Psc based on our orbital solution is 23.10 it 0.24 pc, corresponding 
to an orbital parallax of 43.29 it 0.46 mas, consistent with the Hipparcos result at 3.5% and 
1.7-sigma. 

Based on eclipsing binary systems. Popper (1980) and Andersen (1991) list main-sequence 
linear diameters for stars in the mass range of the 64 Psc components at 1.15 - 1.25 Rq. At 
a distance of approximately 23 pc this corresponds to apparent angular diameters in the range 
of 0.42 - 0.5 mas. This range of apparent diameters is unresolved to our interferometer, with 
expected component visibility moduli (Eq. |^) of approximately 0.99. Consequently, as discussed 
above we have adopted model diameters for the two components based on a fit to Andersen's 
compilation of eclipsing binary data in this mass range. In particular we have constrained our 
solutions to a primary diameter of 0.5 mas (1.24 Rq @ 23 pc) and a secondary diameter 0.45 
(1.12 Rq @ 23 pc) respectively, and included 10% one-sigma errors around these constraints 
into our systematic error computations. We have also ignored any corrections due to stellar limb 
darkening; our data would be highly insensitive to these effects. 

Similarly, because the system is nearly equal mass, hence nearly equal brightness, our V'^ 



data is relatively insensitive to the magnitude difference of the two components (see Hummel et 



al. 1998 for similar comments). Based on our masses and the mass-luminosity models of Henry & 
McCarthy (1992, 1993) one would expect a ii'-magnitude difference of 0.11 and a ^-magnitude 
difference of 0.24. For reference, N79 quoted a best-fit component magnitude difference in 
the visible of 0.16 without a quoted error. An independent fit of our PTI data for a K-band 
magnitude difference between the two yields components 0.141 it 0.076 - well within one sigma 
of the expected i^-band magnitude difference. As noted above for our orbit fit we have chosen 
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to constrain the ii'-magnitude difference at its model value, and included a 5% one-sigma error 
around this constraint in our systematic error computations. However, our fit K-band magnitude 
difference, system distance estimate, and a 64 Psc apparent (CIT) K-band magnitude estimate 
of 3.88 ± 0.07 yields K-band absolute magnitude estimates of 2.75 ± 0.08 and 2.89 ± 0.09 
for the primary and secondary components respectively. These estimates agree well with the 
predictions of the mass-luminosity model of Henry and McCarthy ( Henry &: McCarthy 199^ ), 
which predicts absolute K-band magnitudes of 2.79 and 2.91 for the primary and secondary 
components respectively. 



5. Temporal Variations in Orbital Elements 

N79 argued for the possible presence of an additional, undetected object in the 64 Psc system 
from supposed temporal variations in the systematic velocity, the time of periastron (To), and 
the argument of (primary) periastron (wi). In particular, N79 gives plots of systematic velocity, 
projected Tq, and coi as a function of observation epoch. Assuming the N79 estimate for the 
period of the effects ( 2700 yr), the approximate angular separation between 64 Psc and this third 
component would be ~ 10 arcsec, sufficiently distant to not affect our V'^ measurements on the 
inner pair, and thus not perturb our orbit reconstruction. 

Our interferometric data does not address the 64 Psc systematic velocity, but we can test 
the consistency of the N79 orbital element trends with our orbit solution. Because we note a 
bias in the period estimate could produce the linear effects in projected Tq that N79 claims, the 
comparison is best done with oji. Figure ^ gives a plot of various values of uJi as reported by 
AL76, N79, DM91b, and PTI. Updating the N79 model for ivi, we have fit a linear model to wi 
based on 1918 measurements made at Lick and Mt. Wilson (reduced by N79), AL76 (reanalyzed 
by N79), N79 themselves, and DM91b. As depicted in Figure ^, this model is dominated by the 
estimate determined by N79 from the 1918 Lick/Wilson data, and based on the errors quoted 
by N79 the data will support a broad range of trends for the time-variation of wi - including 
essentially no variation. The updated N79 model for uji evaluated at the epoch of our 1997/1998 
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Physical 


Primary 


Secondary 


Parameter 


Component 


Component 


a (10-2 AU) 


7.371 ± 0.052 (0.052/0.006) 


7.706 ± 0.054 (0.054/0.07) 


Mass (M0) 


1.223 ± 0.021 (0.021/0.002) 


1.170 ± 0.018 (0.018/0.002) 


Sp Type (DM91b) 


F8V 


F8V 


Model Diameter (mas) 


0.5 (const) 


0.45 (const) 


System Distance (pc) 


23.10 ± 0.24 (0.13/0.21) 


Tiorb (mas) 


43.29 ± 0.46 (0.24/0.39) 


Mk (mag) 


2.75 ± 0.08 (0.07/0.04) 


2.89 ± 0.09 (0.08/0.04) 



Table 4: Physical Parameters for 64 Psc. Summarized here are the physical parameters for 
the 64 Psc system as derived from the orbital parameters in Table ^. As for our PTI-derived 
orbital parameters we have quoted both total error and separate contributions from statistical and 
systematic sources (given as (t stat / sys) ■ 
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PTI observations is not in particularly good agreement with our V^-fit and full-fit values for loi of 
202.7 ± 1.9 and 203.56 ± 0.35 respectively (Table 1); it predicts a value for uji of 205.0 ± 2.0 deg. 
Shown in Figure the disagreement (4-sigma) of our full-fit wi determination to this linear model 
central projection would seem to cast doubt on the evolving-tJi hypothesis at a reasonably high 
confidence level. The potential problem with drawing a firm conclusion on the evolution question 
is the inclusion of the older DM91a radial velocities - while the RV data match our recent visibility 
measurements very well, they may well bias the oji estimate away from its current true value. 



6. Companion Search 

Because our 64 Psc orbital studies do not rule out the presence of a companion as N79 suggest, 
we imaged the 64 Psc system with the Near-Infrared Camera (NIRC) on the Keck-I telescope at 
Mauna Kea on 20 October 1998, and again on 28 November 1998. A panel of near-infrared images 
of the 64 Psc system is shown in Figure |^. 

These images were obtained through the J, H and K filters with total integration times of 
60 s - the data from six 10 s exposures being coadded into a memory buffer before being recorded 
to disk. In the J band we took 10 such images while moving the telescope pointing by 10 pixels 
between each exposure to limit the effects of bad pixels on the array. These 10 images were then 
shifted and added. The same procedure was used at H and K band, but only 7 images were 
taken in H and 1 in K. These were sufficient to detect objects down to J = 21.0, H = 20.2 and 
K = 19.0, where the total exposure times are 600 s, 420 s and 60 s for J, H and K respectively. 
Clearly the presence of the bright 64 Psc system in the images makes detecting faint objects 
difficult. By reflecting the images along the central vertical axis (where the star was positioned) 
and subtracting this reflected image from the original, a substantial portion of the star light is 



removed. (See Matthews et al. 1996 for a detailed description of this technique.) This subtracted 
image is used for photometry of the remaining objects. We used the standard SJ9101 (HST 
P525-E, Persson et al. 199^) to calibrate the instrumental magnitudes and determine the flux 



density of the point sources we found in these images. 
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Fig. 3. — Measurements of wi vs. Time for 64 Psc. At least as calculated by N79, measurements 
of uji vs. epoch of observation show a striking linear trend, and N79 suggest the presence of an 
additional, undetected body in the 64 Psc system to explain these effects. The N79 linear model for 
LOi (updated here to include the DM91b determination, but not the PTI result) is not in particularly 
good agreement with the value for wi estimated from our 1997/1998 PTI observations of 64 Psc 
(both separately and in combination with the DM91a radial velocity measurements). Our combined 
interferometry-radial velocity determination of coi is sufficiently accurate to formally exclude the 
N79 evolution hypothesis at high confidence, but it conceivably may be biased by the older radial 
velocities. 
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There are two apparent point sources other than 64 Psc itself in our images, which following 
the WDS conventions we will refer to as 64 Psc D and E (momentarily leaving aside the issue of 
whether these objects are actually assocated with the 64 Psc system). To determine the offsets of 
the D and E objects with respect to 64 Psc we had to determine the position of 64 Psc on the 
detector, which was completely saturated within a radius of 23 pixels of the star. (The pixel scale 
is 0.15".) The star's position was determined by fitting lines through the unsaturated portions 
of the 6 diffraction spikes yielding a position with a conservative accuracy of better than 0.05" 
(a third of a pixel). The first point source, herein 64 Psc D, is in the upper right of the images, 
separated from 64 Psc by 22.7" at a position angle of 250°, measured N through E. This point 
source has the following magnitudes: J = 17.0 it 0.1, H = 16.8 ± 0.1 and K = 16.6 ± 0.1. This 
photometry is consistent with that of a background G7 V star at a distance of roughly 4 kpc. It is 
therefore apparently excluded as a possible companion to 64 Psc. 

The other point source, herein 64 Psc E, is much dimmer than the D companion, and we have 
highlighted it in the figures with a red arrow to indicate its position. 64 Psc E is separated from 
64 Psc by 12.1" at a PA of 60°, and is extremely red; it is undetected in the J image, but visible 
in H and clearly detected at K bands. 64 Psc E shows no obvious signs of extended structure in 
either the H or K- band images. Photometric measurements of 64 Psc E yield the following: J > 
21.0, H = 19.3 ± 0.2, K = 18.5 ± 0.15. If this candidate is physically associated with 64 Psc it 
would have an Mk = 16.7 making it 2.4 magnitudes fainter than Gliese 229B, the coolest known 
brown dwarf and the least luminous condensed object directly detected outside the solar system. 

Because of 64 Psc Es intriguing red color, and the fact that the separation between it and 
64 Psc is so near the calculated separation if the N79 conjecture were correct, we obtained a 
near-infrared spectrum with NIRC on the night of 28 November 1998. 

To obtain the spectrum we used NIRC's 120 line mm~^ grism in conjunction with the 4 pixel 
slit (chosen to match the seeing) and blocking filters which permit only the H and K bands to 
pass through to the detector. The object was positioned on the slit and 14 exposures of 60 s 
integration (1 coadd) were taken, while moving the object along the slit by 10 arcseconds between 
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Fig. 4. — Keck NIRC images of 64 Psc. Left-to- Right, Top-to-Bottom the images are J-band (1.2 
/im), ff-band (1.6 /im), ET-band (2.2 ^m), and a K minus i^-reflected image (see the text for an 
explanation of this last image). The field of view is 32" on a side with a pixel size of 0.15". The 
red arrow indicates the peculiar object, which we call 64 Psc E, whose spectrum is shown in Fig. |^. 
64 Psc Ks J - K color is > 2.5"^. 
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exposures. This permitted us to extract spectra from every exposure. We then observed the 
standard G8III star SAO 92183 to provide a spectral flatfield from the same airmass. SAO 92183 
was smeared across the entire length of the slit to produce a "flatfield" over the whole detector. 
A sky frame was taken identically to be used for subtraction from the reference (SAO 92183) 
spectral image. The data were reduced by summing all the spectral images and dividing by the 
reference spectral image. Then the spectrum of 64 Psc E and that of one of the diffraction spikes 
of 64 Psc were extracted separately. The two were separated by 5" in the spectral images. This 
means that it is unlikely any contamination from 64 Psc is present in the spectrum of 64 Psc E. 
The spectra were extracted by summing the 7 rows over which the signal from each source was 
visible. The extracted spectra were then wavelength calibrated using the well-determined 50% 
transmission points of the HK blocking filter. The flux calibration involved first the removal of 
the Rayleigh-Jeans falloff of the reference G-star used to flat field the data and a matching of the 
i^-band integrated flux to that measured by the photometry. 

The spectra from 1.4 to 2.5 /um of both the candidate companion 64 Psc E and 64 Psc itself 
are shown in Figure ^. The candidate spectrum is extremely red, confirming the photometric 
measurements, and it is substanitally different from the composite spectrum of the two late F-star 
components of 64 Psc. This demonstrates that little or no contamination from 64 Psc itself is 
present in the spectrum of 64 Psc E. While the SNR is low, there are no apparent molecular 
features in the spectrum. In particular the large methane features seen in the spectrum of the 
cool brown dwarf GL 229B (see Oppenheimer et al. 1995, Geballe et al. 1998, pppenheimer et 



al. 1998) are notably absent. 



Based on the colors and proximity of the first (64 Psc D) candidate, and the implied large 
absolute magnitude of the second (64 Psc E) candidate, we think it is unlikely that either object 
is physically associated with the 64 Psc system. However, 64 Psc E is relatively unusual in its 
extreme red colors. Based on initial statistics from the 2MASS survey, L-dwarfs typically have J 



K colors less than 2.1, so that possibility would seem to be unlikely (Kirkpatrick et al. 1999) 



Other possible explanations for the extreme red nature of 64 Psc E is the possibility of a very red 
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Fig. 5. — The near-infrared spectrum of the candidate companion of 64 Psc, 64 Psc E, is shown 
in the top panel. The bottom panel is the composite spectrum of the two F stars that comprise 
64 Psc. That spectrum was taken from the diffraction spike which passed through the slit of the 
spectrograph. The gap in both spectra between 1.9 and 2.0 //m is where the atmospheric absorption 
prevented measurement. We estimate the SNR in the 64 Psc E spectrum is approximately 2 in the 
blue edge of the i?-band and 35 in the red end of the i^-band. 



QSO ( Beichman et al. 199^ ), or a red, compact background galaxy (Becklin et al. 1995, Graham & 
Pey 1991 , [Stanford et al. 19971 ). It is clearly interesting to determine whether the 64 Psc E object 
shows indications of common proper motion with 64 Psc itself, and a higher SNR spectrum of the 
object would be needed to unequivocally establish the nature of 64 Psc E. 



7. Summary and Discussion 

By virtue of our interferometric resolution and the high precision of the DM91a radial 
velocity data we are able to determine accurate masses of the 64 Psc constituents, and an accurate 
system distance. Our data favors the earlier spectroscopic period of N79 over the more recent 
DM91b determination. The approximate determination of relative K-brightness ratio of the two 
nearly-equal brightness 64 Psc constituents (best-fit K-magnitude difference of 0.141 it 0.076) 
prevents us from establishing tight constraints on the absolute magnitudes of the components. 
This is regrettable, but is an inherent limitation of our amplitude-only reconstruction technique 



for nearly equal-brightness systems ( Hummel et al. 1998 ). 



The N79 suggestion of an additional body in the system as an explanation for purported 
orbital parameter variation is intriguing, but our orbital results do not seem to favor this inference. 
However, these considerations suggest high angular resolution, high dynamic range imaging of 
the system to attempt a detection of a companion, and we have imaged the 64 Psc system with 
NIRC at Keck-I. Neither of the two potential companion candidates apparent in our near-infrared 
images seem likely to be physically associated with the 64 Psc system. However, the second of 
these candidates (herein 64 Psc E) would appear to be fairly unusual based on its very red colors 
and spectrum in the near-infrared, and follow-up proper-motion and spectroscopic observations 
would be warranted. 
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